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Q" 1 Abstract 

kt\| In this paper we provide fully covariant proofs of some theorems on shear- 

free perfect fluids. In particular, we explicitly show that any shear-free per- 
fect fluid with the acceleration proportional to the vorticity vector (including 
rS | the simpler case of vanishing acceleration) must be either non-expanding or 

non-rotating. We also show that these results are not necessarily true in the 
Newtonian case, and present an explicit comparison of shear-free dust in New- 
tonian and relativistic theories in order to see where and why the differences 
appear. 

PACS: 04.20.-q, 04.20.Jb, 04.40.Nr, 98.80.Hw 



1 Introduction 

This paper deals with shear-free perfect-fluid solutions of Einstein's field equations. 
The motivation for this study comes, on the one hand, from some studies on kinetic 



theory (see [37| and references therein), and on the other, from the relations between 
relativistic cosmology and Newtonian cosmology Concerning the former, when we 
consider isotropic solutions of the Boltzmann equation, that is, those for which 
there is a timelike congruence with u as the unit tangent vector field such that 
the distribution function has the form f(x a , E) with E = —u a p a (where p a denotes 
the particle momentum), two important results follow: i) The energy-stress tensor 
computed from such a distribution has the perfect-fluid form with respect to u (see 



for instance |33], 35]). ii) The unit tangent vector field is shear- free and in addition 
its expansion 9 and rotation u satisfy u9 = (see the proof in |37j). These results 
led to the formulation of a conjecture whose origin seems to be the Ph. D. thesis 
by Treciokas []36| (see [^7j] for more details). This conjecture can be expressed in 
the following form (here g and p are the energy density and pressure of the perfect 
fluid): 

Conjecture 1: In general relativity, if the velocity vector field of a barotropic 
perfect fluid (g + p ^ and p = p(g)) is shear-free, then either the expansion or the 
rotation of the fluid vanishes. 

While we are still probably a long way from settling the truth or falsity of Conjecture 
1, it is something short of amazing that such a conjecture might be expected at all 
in general relativity. Consider for example the pressure-free (dust) case for which 
Ellis |J5J| showed that o = ==>- Quo = 0. This is a purely local result to which 



no corresponding Newtonian result appears to hold, as counterexamples can be 
explicitly exhibited p2| |. Ellis's theorem holds for arbitrarily weak fields and fluids 



of arbitrarily low density. Why then does the Newtonian approximation fail? 

Knowing whether or not this conjecture is true, or at least to what extent it is 
valid, might be useful in seeking and studying new perfect-fluid solutions of Ein- 
stein's field equations with a shear-free velocity vector field. With respect to this 
subject, there are some interesting studies of shear-free perfect-fluid models to be 



found in g, |9|, [10], [II]. On the other hand, it is important to remark that there are 
many known cases which are shear-free and either rotation-free or expansion-free. 
Some examples are: the Friedmann-Lemaitre-Robertson-Walker space-times (see 



for instance |2lL [23|]), the Godel solution |L9], |21]j , spherically symmetric shear-free 



perfect-fluid solutions with expansion and equation of state |]25|j, Winicour's station- 



ary dust solutions |4(| , and other examples of perfect-fluid stationary solutions with 
barotropic equation of state where the velocity vector field is aligned with a timelike 
Killing vector field. 



The conjecture has been proved in some special cases. As far as we know these 
cases are the following: 



i. Spatially homogeneous space-times: Schucking |3T| studied the case of dust 
(p = 0); Banerji Jl| studied the case with a linear equation of state p = (7 — 1)^ 
with 7 > 1 and 7 7^ 10/9; and finally King and Ellis |M showed the general 
case with g + p > 0. 

ii. Ellis |L5| proved the conjecture for dust (p = 0). 

iii. Treciokas and Ellis [37| showed the case for incoherent radiation (p = l/3g). 



iv. The case in which the acceleration and the vorticity of the perfect fluid are 
parallel was shown by White and Collins 



v. The case in which the magnetic part of the Weyl tensor with respect to the 
velocity vector field vanishes was proved by Collins ||. 

vi. Carminati [|j] showed that for the case of Petrov type N the conjecture also 
holds. 

vii. Lang and Collins |27| proved the same for the case in which the expansion and 
the energy density are functionally dependent {6 = 0(g)). 

viii. Coley M considered the existence of a conformal Killing vector field parallel 
to the velocity vector field u, proving also the conjecture for this case. 

ix. Finally, there are some recent partial results on Petrov type III by Carminati 
and by Carminati and Cyganowski |J. 

It should be stressed, however, that all the above results have been proved using 
either a particular tetrad or coordinate system. In this sense, there is a clear needtl 
for a fully covariant proof of some of the above partial theorems. Such covariant 
proofs are not only desirable on aesthetical grounds, but they may also be useful 
for the deeper understanding of why the theorems hold. They may thereby aid in 
the further development of the subject, perhaps helping in eventually proving or 
disproving the above Conjecture 1. 

In this paper, we present fully covariant proofs of the theorems in two relevant 
cases: when the acceleration of the fluid vanishes [including the above case (ii) 
of dust] and when the acceleration is parallel to the vorticity vector [case (iv)]. 



lr This was emphasized by G.F.R. Ellis in his plenary lecture during the past Indian Inter- 
national Conference on Gravitation and Cosmology (ICGC-95), held in Pune, December 1995. 
Two of us (JMMS and PS) were attending this lecture which aroused our interest on the subject 
independently. 



Analogous covariant proofs can also be given for other cases (see [53]) by using similar 
procedures. In Section 2 we present the main equations and Sections 3 and 4 are 
devoted to the proofs of the two theorems. Finally Section 5 deals with Newtonian 
theory, and by attempting to follow our general relativistic proof, uncovers the 
reason for the failure of the Newtonian limit. In addition, all shear-free Newtonian 
universes are revealed by this analysis. 

2 General results on shear-free perfect fluids. 

Let us consider a perfect fluid with unit velocity vector field u, so that the energy- 
stress tensor readsH 

T ab = OU a U b+P P abi (!) 

where g and p are the energy density and the pressure, respectively, and 

P a b — 9ab + u a u b-. Pab — P{ah) (2) 

is the projector orthogonal to u, which has the standard properties 

P a b P b c = P% PI = 3, P a b u b = 0. (3) 

Let us summarize briefly the main concepts for the study of the kinematics of the 
velocity u. The derivative along u of any tensor quantity with components A b ... b v 
will be denoted by 

u a v a AXX = a%:% . (4) 

Some different names for this derivative are used in the literature: material deriva- 
tive, convective derivative, time derivative with respect to u, etc (see p8| ). Along 
this paper we shall always use the term time propagation along u. A particular 
and interesting case is the time propagation along u of u itself, which is called the 
acceleration vector field of the fluid and is denoted by 

a a = u b V b u a = if. (5) 

^From this definition and taking into account the fact that u is a unit vector field, it 
follows that the acceleration is orthogonal to u and therefore it is a spacelike vector 
field. Moreover, the integral curves of u are geodesic only when its acceleration van- 
ishes, so that the fluid is said to be geodesic if the acceleration vanishes. Physically, 
the acceleration vector field represents the mixed effects of gravitational as well as 
inertial forces (see refs.|3l, 11(3]). 



2 Throughout this paper we use a,b,c,...,h = 1,2,3,4 for spacetime indices, while we use 
i,j, k, . .. = 1,2, 3 for indices in the Newtonian theory, see section 5. 



The spatial part of the covariant derivative of u decomposes in general into its 
irreducible parts with respect to the rotation group (see [|16|, |17[ ^] and references 
therein) as follows 

^b u a + a a w b = ~0P ab + o ah + u ab . (6) 

where 

e = vx, (7) 

a ab = P£PJV d u c - V a „ a ab = a {ab) , a\ = 0, a ab u b = 0, (8) 

W «t = P {a P b]^ d U c , U ab = U [ab] , U ab U b = 0. (9) 

These objects together with the acceleration a form the kinematical quantities of u. 
The commonly used names for them are: expansion rate for 6, shear rate for a ab and 
rotation rate for uo ab . 

We can also introduce the vorticity vector field by 

" a = \n ahcd u b u cd (io) 

from where we see that the vorticity is orthogonal to u, that is, u a u a = 0. Moreover, 
we can invert fllCf ) in order to obtain the rotation rate in terms of the vorticity 

Uab = VabcdU C U d (11) 

where rj abcd is the canonical volume element in the spacetime (see, for instance, 



lj, [H], £25J). ^Prom ([Tj) it follows directly that u ab u b = 0. We also define the 
rotation scalar by 

cu 2 = l -u ah u ab = u a u a , (12) 

and, given that u ab is a spatial tensor (that is, completely orthogonal to u), u 2 is 
non-negative and the following equivalences hold: 

uj = ^^ uj a = ^^ uj ab = 0. (13) 

^From (|Tl"| ) we can immediately get the following standard identities 

u a c u b = uo a u b - uo 2 P b =► < ^ c d ^ = -u> V . (14) 

Finally, let us remark an important fact concerning u ab . Taking into account the 
following expression which comes directly from (|6]) 

V [a « 6] = & [a u b] + uj ba , (15) 

it follows that 

U [c ^a U b] =0 ^^ 0J ab = 0, (16) 



and using Frobenius's theorem (see JT8], |25|), this means that there exist locally 
functions / and h such that u a = hV a f, (or equivalently, u generates orthogonal 
spacelike hypersurfaces), if and only if the rotation vanishes. 
The fluid (or also u) is said to be shear-free if 



On 







(17) 



which we assume from now on. In this case, equation (19) becomes simply 

V fe w a = -0P ab + uj ab - a a u b . 
Then, by using the Ricci identities for u 

( V c V d - V d V c ) u a = R a bcd u b 



;is) 



(19) 



we can compute the time propagation along u of (|18 ), which produces the following 



set of evolution equations [0, [13|, |2(| (we use units with 8irG = c = 1) 



+ ^ 2 -2^ 2 -V a a a + ^(f? + 3p) = 0, 



1 c TD d,'. 



2 
3' 



> c d dx 



-P^-Pfe ^cd + o^ai< ~~ ^[a-ffe] V rf a c 



0. 



£. 



afc 



(20) 

(21) 
(22) 



o^ -^afe + a a a b + E {a E b) ^c a d ~ TT • 

Here _E ab is the so-called electric part of the Weyl tensor with respect to u. This and 
the magnetic part H ab relative to u determine completely the Weyl tensor C acbd , and 
are defined respectively by 



E ab — ^acbd U u 



E ab = E (ab)i E ab U = °> E \ 



0. 



(23) 



H ab = ^C^ul 



H ab = H {ab) , H ab u b = 0, H a a = 0. (24) 



These tensors were first introduced by Matte p9f in the context of the study of 
gravitational radiation. 

Equation ( p0|) is the famous Raychaudhuri equation pOj , and Eq . fl2~l|) can also 
be given in terms of the vorticity as follows 



P a b u b + 2 -6u a - l -v abcd u b V d3 , c 



0. 



(25) 



Eq. (|22| ) comes directly from the time-propagation of the absence of shear. The whole 
set ( pmplj ^2|) contains 9 of the 18 independent components of the Ricci identities 
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I9|). The remaining 9 components are usually called constraint equations and are 



given by [12, 13, 20 



2P ab V b 6 + 3P a b V d u bd + 3^ a fe a 6 = 0, (26) 

V a u a = 2a a w a «=► V abcd u a {V b u cd - ^ cd ) = 0, (27) 

H ab = 2u {aH) - {^ c )P ab + P { °P b fV c u d , (28) 



where in order to write down the magnetic part H ab in the form (|2~8|) we have used 
the identity 

u ab V c u bc = -u\ + (a 6 ^R + P a b Lu c V c cu b - P a b uj c V b u c . (29) 

Finally, the conservation equations V b T ab = for the energy-stress tensor ([!]) 
read 

q + (q + p)6 = 0, (30) 

(g + p)a a + P ab V bP = 0. (31) 

Now, we are ready to present the fully covariant proofs of the main theorems. 

— » 

3 Case without acceleration: a = 0. 

In this section we present the coordinate- and tetrad-free proof of Conjecture [1] for 
the case of vanishing acceleration 

a = 0, (32) 

which includes the particular case of dust because of (pl|). First of all, we present 
some lemmas which will help us in the proof of the theorem. 

Lemma 1: If there exists a function f satisfying P ab V b f = then either f =const. 
or the rotation vanishes. 

Proof: There are two possibilities, either / vanishes or not. If it does, then 
/ =const. clearly. If it does not, then we can write 

u a = -jVJ. (33) 

But then, as we have seen in the comments following equation flT6"|), this implies 
u ab = 0, proving the lemma. 

Lemma 2: If the perfect fluid is geodesic then either the pressure p is constant or 
the rotation vanishes. 



Proof: The proof of this lemma comes from the energy-stress conservation equa- 
tion ([SID , which by using fl32| ) becomes P ab V b p = 0. Then Lemma 1 implies the 
result. 

Lemma 3: If the perfect fluid is geodesic and shear-free, and there exist constants 
C\ and c 2 with c 2 ^ such that 

g={ Cl -l)p + c 2 uj 2 , (34) 

then either the rotation or the expansion vanishes. 

Proof: First of all, the geodesic condition implies that the time-propagation 



equations (|21|) and ( P5|) for the rotation and the vorticity reduce to 



P a c P b d "c d = " ab = ~8"a b , (35) 

p a b u b = u a = -~9u a => uj = --duo. (36) 

o o 

^From Lemma 2 either the rotation vanishes or the pressure is constant. Thus, we 
need only consider the case p =const. The time propagation of ( |34"|) then gives, on 
using (J30|), and (|36|) , 

(c ip _ ^ = o. (37) 

If 9 = we are done. If the term in brackets vanishes, then uj is constant, and (|36|) 
gives 9 uj = as required. We now pass to the proof of the theorem. 



Theorem 1: Every shear-free and geodesic perfect fluid must be either expansion- 
free or rotation- free, i.e. 

a ab = 0, a b = => uo9 = 0. 



This theorem was first shown by Ellis [0 for the dust case (p = 0) and it was 



completed for the case of constant pressure by White and Collins |39(|. Both theorems 



were proved using the so-called tetrad formalism (see, for instance, p8 



Proof: From Lemma 2 it is clearly only necessary to consider the case p=const. 



Moreover, from ( fLq) restricted to the case with (|32D we get immediately, for any 
function /, the following identity 

(P ab V b f) ■ = P ab VJ + uj ab VJ - 6 -P ab VJ . (38) 



Using this identity together with the relations ( |35| - |36| ), the time propagation of 
Equation (|26f) restricted to the case a b = results in 



b, ,c\ 



2P a b V b g - 13P> c V b ^ c - 3P a b u c V c u b = 



(39) 



On using (p(f) contracted with u ac and the identity (|29|) , this equation can be rewrit- 
ten as 

Py b g - 8coPy b cu + uj b a V b 9 = 0. (40) 

This key equation provides an algebraic relation between the spatial gradients of 
density, rotation and expansion. By contracting it with u> ca we also get 



uy b6 - 8uuy b u + u a c uy b 9 = o. 



(41) 



The next step is the computation of the time propagation of Equation ([K]). Making 
use of ( J35l - |3~6D , (|38"D, (|30"D, (^TJ) and the Raychaudhuri equation (|20|) this leads to 



q + p - y ^ 2 ) p a b v b e = \o^y b e + e -py bB . 



(42) 



Time-propagating again this relation and taking into account the key equations fl4H| ) 
and (f|l]) together with (|35[ - |36"1) , ([38]), (p0[), and ([42]) itself we arrive at 



i^-i> a v+» 



112 2 5. 



^ 6 v^ 



" ( 1^ - r 2 



H^) ^^ + l 0UJ >c V » 6 ~ \"a^y b e = . (43) 



Multiplying this by 9 in order to use ([l"2) again, and re-ordering terms it follows 
that 



1(q + P- y- 2 ) ^V b 9 - y^y b 9 + (y + - 2 " ^) "^ h 9 



+ 



32 
~9~ 



flV+^ + p-^V) 



^±^-2^ 






P a b V b 9 = . (44) 



Contracting with uo a we find the simple relation 



lev+^+p-H^V^-^ 



o/V 6 # = 0, 



(45) 



so that two different possibilities appear: 
1) If uj b V b 9 ^ 0, then 



feV+( e +p-f</)(^-2^ 



0. 



(46) 



The time propagation of this equation leads on using Eqns. (|20|), fl30|), fl36|) and (|46| ) 
itself, either to = (in which case obviously 9u> = 0) or 

—u 4 - 2qlu 2 - ^pco 2 +p(g + p) = 0. (47) 

Once more time-propagating this relation we finally arrive at 

7 40 2 \ 2 



e +-p-— u 2 ) u/ = 0. (48) 

Thus, either u = and we are done, or the term in brackets vanishes in which case 
uj9 = by Lemma 3. 

2) If u b V b 6 = 0, then Equation @ together with the identities (|T|) lead to 



6 (29cu 2 -Q(g + p))uJ a b V b e 



5 V .2 , fn.2 Q + 3P\ / on .2 



2 u; 2 + ( 2c^ 2 - ^^ ) (29^ 2 - 6( e + p) 



P a "V^ . (49) 



But the two vectors P a fe V fe ^ and u; a 6 V 6 # are orthogonal to each other, as follows 
immediately from (|i~l|) , 

w 6 V 6 ^ = ^cP c6 V 6 = r lttea jJ i u'P d V# . (50) 

Hence the left-hand side of relation fl49"D must vanish, and there are three possibil- 
ities. Either (a) 9 = 0, in which case the theorem is proved, or (b) the term in 
brackets vanishes whence the theorem follows from Lemma 3, or (c) uj b V b 8 = 0. In 
this case using the condition defining case (2), namely u b V b 8 = 0, it follows that 
either P b ^ b 9 = or u b = 0. If the latter holds then u = and the theorem is 
proved, while the former implies through (]42|) and (|40D that P a 6 V b £ = P a 6 V fe co> = 0. 
Then setting / = uj in Lemma 1 gives u is constant and Equation ( |56| ) implies that 
0a; = 0. Hence the theorem holds in this case too. 

This finishes the proof of Theorem |l| 

It must be remarked that Theorem 1 does not need the assumptions of a ba- 
rotropic equation of state and g + p ^ 0. In particular, this theorem also holds 
for Einstein's spaces (including vacuum) in which there is a geodesic shear-free unit 
timelike vector field. 

Corollary 1.1: If u a is a shear- free geodesic vector field in any Einstein space 
(Gab = ^9ab), then 9u = 0. 

Proof: Essentially this is the case 

q = —p = —A = const. 

10 



Equation (|45| ) is replaced by the simpler 

"2 



^ 



-e 2 + P + 2u 2 



Lu b V b 6 = 0. (51) 



If the expression in brackets vanishes we have by essentially the same discussion as 
given in case (1) above, 

g + 2cu 2 = 

whence 9uj = follows by Lemma 3. On the other hand if u; b Vb# = then the 
proof follows exactly as for case (2), with appropriate simplifications for example in 
Equation 



4 Case in which the acceleration and the vortic- 
ity vector fields are parallel: a || lj. 

In this section we study the case when the acceleration vector field a and the vorticity 
vector field u are parallel, that is to say, when 

E = ipi3, ^^0. (52) 

Taking into account Theorem 1 we will consider ip as an arbitrary non-vanishing 
function. This also means by ( |3~T1) that, whenever there is a barotropic equation of 
state p = p(g), we have 

p 7^ 0, where p = dp(g)/dg. (53) 

As in the previous section, we collect some useful lemmas before we prove the the- 
orem. 

Lemma 4: For every shear-free perfect fluid with a barotropic equation of state 
and g + p ^ 0, we have 

P [6 c P a fV c a d = p'6u ah) (54) 



P a b a b 



v - q - {q + p)— 

3 p 



6a a + uj a b a b + p'P ab V b e. (55) 



Proof: If we introduce the equation of state p = p(g) into relation ([H|) and use 
fl3"0p we obtain 



a a = p' 0u a -(g + p)~ l V a p = p' (0u a - (g + p)~ l V a g) , (56) 

from which it is easily derived that 

P [b c PJl V c a d = p' 6 P [6 c P a f V c u d , (57) 

and now using (fOJ) we arrive at (^). On the other hand, the time-propagation of 
fl56"D , on using ( p6|) and (|3TD, leads directly to (|55|). 
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Lemma 5: For every shear-free perfect fluid with a barotropic equation of state 
and q + p 7^ , the time propagation of the rotation and vorticity are 



PaPb u cd 



"•■'' l p'--)9co a 



P 



9uj„ 



P> 



u> 



P 



! ]e «- 



(58) 
(59) 



Proof: The proof is straightforward on introducing 
fl25|), respectively. 



into equations (|2T| ) and 



Theorem 2: Every shear- free perfect fluid with a barotropic equation of state and 
Q + P 7^ i> n which the acceleration and vorticity vector fields are parallel is either 
expansion-free or rotation-free. 



The proof of this theorem was given by White and Collins |39| by using the tetrad 
formalism |2"8" |. Here, we present a very simple covariant proof of this result. 

Proof: ^From the expression flTT| ) of the rotation rate u ab in terms of ci;, taking 
into account ([52"D , and using equation fl54|) of lemma 4 we can write 



P ab V d u; b 



VabcdU^ 



-v 



-r 1 (^v k f + 2 P 'eu a ) , 



(60) 



which allows us to express Equation ( p6[) in the form 



p ab v b e = ^- l ( -u ab v b *p + p'Ou 1 



(61) 



Introducing this into expression (|55|) of Lemma 4 and using tu^a 6 = ip uj a b uj b = we 
obtain 



a • b 



Plet 



"■-l-^ + 4 +3 \i 



-1„', .ah* 



9& a + -ij-'p'u a0 V b 2P 



(62) 



However another expression for the time propagation of the acceleration vector field 
a can be obtained by putting the assumption (^) into Equation ( [5*9]) of lemma 5. 
The result is 



P%k h 



P 



2 -)o + t 



a a . 



,ab\ 



Comparison of Equations (|62]) and (|63"D leads to 



■0 



^ a0 V^ = 0, 



6ip 



(63) 

(64) 
(65) 
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and the first of these, through Equation (pi), leads in turn to the following simple 
formula 

V a 9 = -9u a + 3iJ>- 1 i/eu a . (66) 

Another consequence of (B3j) is that relation (BDJ) simplifies to 



P*V#> 



bd 



-2 P -9u a 



(67) 



and on contracting with uj^ and using the identity (E9|) we get 



P°»u c V [c u 



b] 



0. 



(68) 



Considering this expression together with the magnetic part of the Weyl tensor (T2S|) 
and the Ricci identities for u ab , the time propagation of fl6T|) gives 



e = ?>( p -\ e 2 



>. 



(69) 



Putting this into equation fl66|) and applying the integrability conditions V^Vyfl = 
for 9, we obtain that either 9 = 0, which proves the theorem, or 



p_ 

V> 2 



u a VJ 



*&^< 



LU 



(70) 



On the other hand, introducing (p9|) into the Raychaudhuri equation ( P0|) and using 
the assumptions of this case we obtain another expression for u; a V a ^, 



^ a v a v 



■KM 



5)2 _ 2 (1 + ^2)^2 + 



£> + 3p 



(71) 



Comparing these two equations we arrive at the relation 

2 



V> 2 



n!) +^ 



^-2(l + ^ 2 )^ + ^±^ P 



3 ? 



>, 



1 r J. y 



a; 



(72) 



whose time propagation provides an expression containing p'" . On the other hand, 
the ui a V a and ui ab V h derivatives of ([F^) lead to two new equations containing also 
p'" . By eliminating p'" from these three equations two different cases appear 
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Case i): 



Case ii): 



2(1+V>V + 



^u a V a u 2 = 2 



3 P - 



1 / y 

%p I 6 p 



i>- 



v'-^W 



u ab V h u 2 



0. 



(73) 



(74) 



In case i), substituting ( [73] ) into ( [72] ) we obtain 

2 



Hi) ^ 



o 2 



Q + 3p 



and time-propagation of this equation leads to 

e(g + p)(l + 3p') = 0. 



(75) 



(76) 



Here, the second factor cannot vanish because of the assumptions of the theorem. 
If the third factor vanished, Eqns. (|7^) and ([75|) would imply 1 + ip 2 = 0, which is 
impossible. Therefore, in case i) we must have 8 = 0, proving the theorem for this 
case. 

In case ii) the time propagation of Eqns. fl74j) does not give additional informa- 
tion. Consider instead the following Ricci identities for the vorticity 



2V [a V 6] u; c = R 



dab' 



U> 



v c (u% 



u 



V c u b V b u c -u b V b V c u c = R bd uj b u d , (77) 



Taking into account that the combination of (68) and ( [74] ) gives 



u b V b u a = uo 2 



r" + ( p '-l^' 



(78) 



and considering the relations obtained until now, Equation ([77]) becomes 



cd\ 



e + p )j* + w ab W ab = 0, W ab = Q ( c a Q b y c uj d - -Q ab Q cd V c u d 



(aVb) 



(79) 



where Q ab = uj 2 uJ a 2 uJ bc is the projector orthogonal to both u and ou. Taking into 
account that 

Q ab = 0, (80) 

the time propagation of ( [79]) leads again to the simple relation (|76|) . If 6 = we are 
done. Finally, if p' = —1/3, then the time propagation of the relation ([72]) leads, 
since 1 + ip 2 ^ 0, to 9uj = 0, which finishes the proof. 
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5 Comparison with the Newtonian case. 

Newtonian cosmology 

By a Newtonian cosmology is meant 

1. A manifold M = E 3 x R, where E 3 is Euclidean 3-space. For every event (x, t) 
in M the fourth coordinate t registers its absolute time coordinate. 

2. An assignment of three functions g{x,t), p(x, t) and 0(x, t) on M, called re- 
spectively density, pressure and gravitational potential. 

3. A differentiable 3- vector field v(x, t) (the velocity vector field) , defined on each 
Euclidean space t = const, whose components depend differentiably on t. 

4. The following equations hold relating the quantities defined above (summation 
convention adopted even though all indices are subscripts, and units have been 
changed in order to restore G) : 

. _dg _dg 

Q= -n = ^I + V iQ,i = ~Q V i,i ( 81 ) 

at at 

dvi dvi . 1 . N 

Vi = -TT = -TTT + VjVij = -<p yi p^ (82) 

at at g 

0,« = AnGg. (83) 

Homogeneous Newtonian cosmologies 



A Newtonian cosmology will be called homogeneous |34[] if g and p have no spatial 
dependence, g = g(t), p = p(t) and the velocity vector field depends linearly on the 
spatial coordinates, 

Vi = Vijtyxj (84) 

for some matrix of components [Vij] depending only on time. From (|82]) it follows 
that <j) t i also is linearly dependent on Xi, 

<f>,i = -fij(t)xj (85) 

where 

fij = % + V ik V kj . (86) 

Decomposing uy into its standard irreducible parts 



1 
! 3' 



V i,j = 0Q$ij + a ij + W ij ( 87 ) 
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where 



j i,i 



expansion 



aij = -(v iyj + v jti ) - -95ij = shear (89) 

LOij = -(vij - v jyi ) = rotation, (90) 
it is found in the homogeneous case that 

= V k k = 0(t), (Tij = V(ij) - -OSij = <Tij(t), LOij = V[ij] = U)ij(t). (91) 



Uij = —zOvij ( 93 ) 



Now suppose that one sets a^ = in the homogeneous case, then Eqns (pl|-|83|) 
reduce to 

q = -0q (92) 

2 

3 

9 = --9 2 + 2u 2 - 4nGp (94) 

2 
fij = u) ik u kj + -(u 2 - 2nGg)5 ij (95) 

where u 2 = |u>yUfy-. 

Define a variable R = R(t) such that 9 = 3R/R, then it is a straightforward 
matter to integrate (|9^) and (|93| ) 

q = g(t ) R~\ oj i:j = u i:j (t ) R- 2 (96) 



where R(to) = 1 and substitution in (j94|) gives 



5 -4ttG g(to) 2oo 2 (t ) 

R = -3~"^" + "3^~- (9?) 

This integrates to give the Heckmann-Schiicking equation |^2| 

3/? 3i? 2 V ; 

where C is an arbitrary constant. Solutions of this equation represent shear-free 
Newtonian cosmologies which are in general both expanding (R ^ 0) and rotating 
(a; 2 (t ) 7^ 0). The gravitational potential is simply read off from (|95l) using fl85"l) . A 
point of interest is that whenever w 2 (£o) 7^ there is no singularity since R(t) always 
has a minimum value Rq > where R = 0. Thus the big-bang singularity of the 
Newtonian equivalent to the standard FLRW models of general relativity (namely, 
the case u 2 = 0), is easily avoided by giving the model an arbitrarily small amount 
of rotation. This is quite contrary to the case in general relativity where singularities 



are generic to all spatially homogeneous models J2T 
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Newtonian version of our proof 

These homogeneous shear-free solutions are independent of the pressure, which may 
as well be set equal to zero (or constant). The question naturally arises as to why 
there are shear-free dust solutions in Newtonian theory having 9uj ^ 0, when Ellis's 
theorem ensures that none exist in general relativity This difference in the two 
theories is both surprising and interesting, since Ellis's is a purely local result and 
is completely independent of the strength of the gravitational field. One would then 
expect it to hold in the weak field Newtonian limit, yet clearly it does not. What 
then is going wrong in this limit? 

To investigate this question we propose to follow the proof of Ellis's theorem 
given here in Section 2, but try to apply it as closely as possible to Newtonian 
cosmology in an attempt to pinpoint exactly where it is that the proof fails. 

Consider a general pressure-free (p = 0) Newtonian cosmology (not necessarily 
homogeneous) having a^ = 0, i.e. 

Vi = ~<f>,i (100) 

where the Newtonian time evolution operator is given by 



d d d 

dt dt l dj 



and 

<f>,u = \q- (101) 



In this last equation, fl83|) has been recast in units such that 8wG = 1 which has 
been done in order to bring the Newtonian and general relativistic equations into 
closer comparison. Eqns. ( PU| ) and ( |100[ ) correspond to the earlier equations ( P~S| ) 
and (0), while Eqn. (|S1~D reads 

q+q6 = 0, (102) 

which is identical with the relativistic equation ( J30|) when p — 0. 
A useful identity is 

d d d d 1 d d 

dt dxj dxj dt 3 dx; %3 dxj 



which takes the place of (|38|). Now perform the time evolution of Equation (|99"D and 

use ( |103|) to arrive at 

2 
Vij = —^Quij, (104) 
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and 



E 



v 



■j j 



\ S -\?^, 



1 2 1 

J. r 2 c 2 c 

7;(f>,kkOij = LUik^jk — 77^ Oij = —LUiLUj + -U 0,, 



'.I- 



(105) 



(106) 



where u>i = htijk^jk- These equations are clearly Newtonian versions of the rela- 



tivistic equations (|21|) 
equivalent of (BBT), i.e. 



and (p2|) when a b = 0. ^From ( |104|) there follows the 
2 



u 



-ecu. 



Applying u 



i,jk 



Vi,kj to (|99|) gives the cyclic identity 



and 



Ujk 



jk,i 



-(^,fe^i 



,•& 



,j<>ikj 



Contracting Eq. ( |108| ) with Uij gives 

LOUJh 



^ij^ikj 



(107) 

(108) 
(109) 

(110) 



while contraction of ( |109[ ) over r/' gives the Newtonian equivalent of the important 

equation (p6j) , 

3 

0,* = A<' (HI) 



Finally contract Eq. ( |109| ) with ufy- and use ( |110|) to give 

3a; u i . 



ujijOj 



(112) 



Now perform the time evolution of ( |111| ), apply ( |103j ) to both sides and use the 
identities (gOg) , ([T05D , ( |Tl~0D , (|TJ) and (|TT2D , results in the key equation 



Q,k-8w<x) )k + Uki9 t i = Q (113) 

which is basically the same as Equation (pE0|). Actually the Newtonian case is rather 
stronger than that in general relativity since the further identity ( |1 1 2j) , which has 
no relativistic equivalent, results in 



g t k = 5uu tl 



-^ki9,i- 



(114) 



At first sight this is a rather surprising result, since Ellis's theorem is a theorem 
concerning restrictions on solutions of general relativity, so one might expect that 
Newtonian theory should place less stringent restrictions on its solutions. 

In any case one can continue in exactly the same way as the proof given for 
Theorem 1 to arrive at a Newtonian version of Equation (f0|), and its contracted 
version 

0. (115) 



32 /)2 2 , ( 16 2\ (Q o 2 

■9 uj + I o — —uj I I — — 2a; 
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Wi9j 



Newtonian equivalent of Lemmas 1 and 3 

In the final stages of the proof of Theorem 1 use is made of Lemmas 1 and 3. The 
Newtonian equivalent of Lemma 3, with p = is 

Lemma 3': If g = kuj 2 for a non-zero constant k, then 9u = 0. 

This lemma holds in Newtonian cosmology, for on using ( |102[ ) and ( |107| ) 

4 
= —do = 2ku>u = — kOuj 2 . 



3 



Hence 



\k6uj 2 = 



which proves the result. 

Lemma 1 however is a different matter. The obvious Newtonian translation 
would be the following: 

If a function /(x, t) has zero gradient, f^ = then f is constant or uoij = 0. 

Of course a function whose gradient vanishes may still be dependent on time, / = 
/(£). This can have no bearing whatsoever on the rotation of the velocity vector 
field. Hence as it stands, this statement cannot be true. It is worth trying to 
understand why the limiting process does not apply for this lemma. 

Suppose we have a metric g a b = Vab + h a b having energy-stress tensor T a b given 
by Equation ([I]) where (restoring c) 

u a =^,-l + o[^y \h ab \=0^pj, (116) 

where v 2 = v(Vi. The Newtonian limit is the limit of all equations found on letting 
c — > oo. For example, the rotation is given by 



u>. 



i.i 



]^ cu \iJ] = v \iJb ( 117 ) 



as to be expected. However the equation Pj'Vbf = reduces in its first three 
components to 

f,i = ~-Jvi (H8) 

c z 

■ df 
where / = — — V f (Vi. The limit as c — > oo of this equation is clearly f i = from 

at ' 

which no conclusions can be made concerning u>ij = vuj-i. In general relativity a 

conclusion can be reached because Equation ([33|) couples u a and V a /, but in taking 

the Newtonian limit v^ and f }i become decoupled. 

Physically one understands it like this. If uj a i, = then the 4-velocity field u a 

is hypersurface orthogonal, and in a sense there is a universal time coordinate t 
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defined such that u a oc V a t (in the lemma, / could act as such a time coordinate). 
However in Newtonian theory such a universal time coordinate is always defined 
and the condition uj^ = is equivalent to Vi = 0^ for some scalar field <fi having 
in general nothing whatsoever to do with the universal time. For example, <fi may 
just depend on the spatial coordinates x, y, z. Lemma 1 is truly a relativistic result, 
having no Newtonian equivalent whatsoever. 

Final stages of the proof 

Returning to the proof of Theorem 1, we attempt to follow the argument given in 
the final stages. Consider the two cases arising from Equation ( |1 1 5| ) . 

Case 1'): io$j, ^ 0. In place of Equation (f47|) the corresponding Newtonian 
argument arrives at 

— u 2 - q\ uj 2 = 0. (119) 

Lemma 3' immediately applies, showing that 9u = 

Case 2') lOiQ^ = 0. The argument can again be continued as in Theorem 1. The 
Newtonian version of (ff9|) is 



-{2$u 2 -Q Q )uj i3 6j 



|W+(W-|)(29^ 2 -6£ 



(120) 



Now clearly uJijOj is orthogonal to 9 t i (since UijO^O^ = and u%j = —uiji), whence 
the left hand side of ( |120[ ) vanishes. If 9 = the proof is over, and by Lemma 3' 
this is also true if 29cj 2 — 6g = 0. If 0Oij9j = then eijk^j9,k — 0, and combined 
with uOiO^i = it follows that Ui = or 9 t i = 0. The former implies u = and we 
are done, while the latter implies by ( |114| ) that 



e,i = w,i = 0. (121) 

However it is not possible to proceed any further. It is at this very last step that 
our attempt at a proof comes to an end since there is no equivalent to Lemma 1. 
The limiting process has broken down for the reasons outlined above. 

All Newtonian shear-free universes 

To conclude this Newtonian discussion, we propose to show that all pressure- free, 
shear-free Newtonian universes are in fact homogeneous. If a = then the above 
discussion has shown that 9u = or 9 t i = 0. In either case we see from Equation 
( gig) that 

Q,i = 0,i = u ti = 0, (122) 

whence 

g = g(t), 9 = 9{t), w = u(t). (123) 
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^From (^), using ( |122| ) and Vijk = Vi,kj we see that u>ij t k = t^ikj, an d substituting 
in Equation (|108|) gives uij^i = 0, i.e. tOjk = ojjkif). Hence 

v i,j = ^(t) S ij + Uijit) = Vij(t) 

from which it follows that 

Vi = VijXj where Vij = -9(t)8ij + Uij(t) = Vij(t). 

All conditions for a homogeneous cosmology are now satisfied, and it follows that 
all shear-free Newtonian cosmologies are homogeneous (and are therefore in the 
Heckmann-Schiicking class). What this means is that although Ellis's theorem 
does not hold in Newtonian cosmology, it "nearly" holds in the sense that the only 
solutions which violate the theorem are homogeneous. 



6 Conclusions 

Ellis's theorem and its various extensions have up till now always required the use of 
tetrads in their proof. In this paper we have shown that it is possible to prove such 
results by methods which are totally covariant in character, and which do not hinge 
on setting up specific tetrads or coordinate systems. While in this paper we do not 
go further than showing already known results, our proofs of these theorems have 
considerable elegance and lend themselves readily for extension to more difficult 
cases. 

Another advantage of our proofs has been the light they shed on the puzzle of 
the Newtonian equivalence. It seems that the basic reason these results do not hold 
in Newtonian cosmology is that there is already a well-defined universal time in that 
theory. In relativistic cosmologies, there will only be a universal time, defined for 
example by the energy-stress tensor, if the rotation of the 4-velocity field vanishes. 
Thus a rotating universe can have no non-constant scalar field whose spatial gradient 
with respect to the 4-velocity vanishes (the content of Lemma 1). However in the 
Newtonian limit this statement is quite untrue, since there exist functions of the 
universal time alone, even when the 3-velocity is rotating. 



21 



References 

[1] Banerji S 1968 Prog. Theor. Phys. 39 365. 

[2] Barnes A 1973 GRG. 4 105. 

[3] Bel L 1989 Recent developments in gravitation. Proceedings of the Relativ- 
ity meeting- 1989 Edited by E Verdaguer, J Garriga and J Cespedes (World 
Scientific, Singapore) 47. 

[4] Carminati J 1988 GRG. 20 1239. 

[5] Carminati J 1990 J. Math. Phys. 31 2434. 

[6] Carminati J and Cyganowski S 1996 Class. Quantum Grav. 13 1805. 

[7] Coley A A 1991 Class. Quantum Grav. 8 955. 
Collins C B 1984 J. Math. Phys. 25 995. 

[9] Collins C B 1988 GRG. 20 847. 

[10] Collins C B and Wainwright J 1983 Phys. Rev. D 27 1209. 

[11] Collins C B and White A J 1984 J. Math. Phys. 25 1460. 

[12] Ehlers J 1961 Abh. Akad. Wiss. und Lit. Mainz. Nat. Kl. 11 792. 

[13] Ehlers J 1993 GRG. 25. 1225 (English translation of @). 



[14] Eisenhart L P 1949 Riemannian geometry (Princeton University Press, 
Princeton). 

[15] Ellis G F R 1967 J. Math. Phys. 8 1171. 

[16] Ellis G F R 1971 General relativity and cosmology. Proceedings of the inter- 
national school of physics "Enrico Fermi" Edited by R K Sachs (Academic 
Press, New York) 104. 

[17] Ellis G F R 1973 Cargese lectures in physics Edited by E Schatzman (Gordon 
& Breach, New York) 1. 

[18] Flanders H 1963 Differential forms with applications to the physical sciences 
(Dover publications, New York). 

[19] GSdel K 1949. Rev. Mod. Phys. 21 447. 

[20] Hawking S W 1966. Ap. J. 145 544. 

22 



[21] Hawking S W and Ellis G F R 1973 The large scale structure of space-time 
(Cambrigde University Press, Cambridge). 

[22] Heckmann O and Schiicking E 1959 Handbuch der Physik LIII Edited by S. 
Fliigge. (Springer Verlag, Berlin-Gottingen-Heidelberg) 489. 

[23] Islam J N 1992 An introduction to mathematical cosmology (Cambridge Uni- 
versity Press, Cambridge). 

[24] King A R and Ellis G F R 1973 Commun. Math. Phys. 31 209. 

[25] Kramer D, Stephani H, Herlt E and MacCallum M A H 1980 Exact solutions 
of Einstein's field equations (Cambrigde University Press, Cambridge). 

[26] Kristian J and Sachs R K 1966 Ap. J. 143 379. 

[27] Lang J M and Collins C B 1988 GRG. 20 683. 

[28] MacCallum M A H 1973 Cargese lectures in physics Edited by E Schatzman 
(Gordon & Breach, New York) 61. 

[29] Matte A 1953 Canadian. Jl. Math. 5 1. 

[30] Raychaudhuri A 1955 Phys. Rev. 98 1123. 

[31] Schiicking, E 1957 Naturwiss. 19 507. 

[32] Sopuerta C F 1996 Ph. D. Thesis (Universitat de Barcelona). 

[33] Synge J L 1957 The relativistic gas (North-Holland Publishing, Amsterdam). 

[34] Szekeres P and Rankin R 1977 Australian Math. Soc. B20 114. 

[35] Tauber G E and Weinberg J W 1961 Phys. Rev. 122 1342. 

[36] Treciokas R 1972 Ph. D. Thesis (University of Cambridge). 

[37] Treciokas R and Ellis G F R 1971 Commun. Math. Phys. 23 1. 

[38] Truesdell C A and Toupin R 1960 Handbuch der Physik III /I Edited by S. 
Fliigge (Springer Verlag, Berlin-Gottingen-Heidelberg) 226. 

[39] White A J and Collins C B 1984 J. Math. Phys. 25 332. 

[40] Winicour J 1975 J. Math. Phys. 16 1806. 



23 



